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Abstract 

The first step in the generalization of the classical theory of homogeneous equations 
to the case of arbitrary support is to consider algebraic systems with multihomo- 
geneous structure. We propose constructive methods for resultant matrices in the 
entire spectrum of resultant formulae, ranging from pure Sylvester to pure Bezout 
types, and including matrices of hybrid type of these two. Our approach makes heavy 
use of the combinatorics of multihomogeneous systems, inspired by and generaliz- 
ing certain joint results by Zelevinsky, and Sturmfels or Weyman [SZ94], [WZ94]. 
One contribution is to provide conditions and algorithmic tools so as to classify 
and construct the smallest possible determinantal formulae for multihomogeneous 
resultants. Whenever such formulae exist, we specify the underlying complexes so as 
to make the resultant matrix explicit. We also examine the smallest Sylvester-type 
matrices, generically of full rank, which yield a multiple of the resultant. The last 
contribution is to characterize the systems that admit a purely Bezout-type matrix 
and show a bijection of such matrices with the permutations of the variable groups. 
Interestingly, it is the same class of systems admitting an optimal Sylvester-type 
formula. We conclude with examples showing the kinds of matrices that may be 
encountered, and illustrations of our Maple implementation. 
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1 Introduction 



Resultants provide efficient ways for studying and solving polynomial sys- 
tems by means of their matrices. This paper considers the sparse (or toric) 
resultant, which exploits a priori knowledge on the support of the equations. 
We concentrate on unmixed (i.e. with identical supports) systems where the 
variables can be partitioned into groups so that every polynomial is homoge- 
neous in each group. Such polynomials, and the resulting systems, are called 
multihomogeneous. Multihomogeneous structure is a first step away from the 
classical theory of homogeneous systems towards fully exploiting arbitrary 
sparse structure. Multihomogeneous systems are encountered in several areas 
including geometric modeling (e.g. [CGZ98], [Sax97], [ZhaOO]), game theory 
and computational economics [MM97] . 

Known sparse resultant matrices are of different types. On the one end of the 
spectrum are the pure Sylvester-type matrices, where the polynomial coeffi- 
cients fill in the nonzero entries of the matrix; such is the coefficient matrix of 
linear systems, Sylvester's matrix for univariate polynomials, and Macaulay's 
matrix for homogeneous systems. On the other end are the pure Bezout-type 
matrices, i.e. matrices where the coefficients of the Bezoutian associated to 
the input polynomials fill in the nonzero entries of the matrix, whereas hy- 
brid matrices, such as Dixon's, contain blocks of both pure types. The ex- 
amples in Section 7.2 show the intricacy of such matrices. Hence the interest 
to describe them in advance in terms of combinatorial data, which allows for 
a structured matrix representation, based on quasi-Toeplitz or quasi-Hankel 
structure [EP02],[MP00]. 

Our work builds on [WZ94] and their study of multihomogeneous systems 
through the determinant of a resultant complex. First, we give precise degree 
vectors together with algorithmic methods for identifying and constructing 
determinantal formulae for the sparse resultant, i.e. matrices whose determi- 
nant equals the sparse resultant. The underlying resultant complex is made 
explicit and computational tools are derived in order to produce the small- 
est such formula. Second, we describe and construct the smallest possible 
pure Sylvester matrices, thus generalizing the results of [SZ94] and [GKZ94, 
Sec. 13.2, Prop. 2.2], already present in the interesting paper [McC33] pointed 
by one of the referees. The corresponding systems include all systems for which 
exact Sylvester- type matrices are known. We consider more general Sylvester- 
type matrices, and show that in the search of small formulae, these more 
general matrices are not crucial. The third contribution of this paper is to 
offer sufficient and necessary conditions for systems to admit purely Bezout 
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determinantal formulae, thus generalizing a rcsTilt from [CKOO]. It turns out 
that these are precisely the same systems admitting optimal Sylvestcr-type 
formulae, and this is nothing but a special case of complexes with only two 
non- vanishing cohomologies. We also show a bijection of such matrices with 
the permutations on {1, . . . , r}, where r stands for the number of the variable 
groups. While constructing explicit Bczout-type formulae, we derive a precise 
description of the support of the Bezoutian polynomial. 

The complex with terms K^{m) described in the next section is known as 
the Weyman complex. For any choice of dimensions, of degrees of the input 
equations and of an integer vector m, the multihomogeneous resultant equals 
the determinant of the Weyman complex (for the corresponding monomial 
basis at each of the terms) , which can be expressed as a quotient of products 
of subdeterminants extracted from the differentials in the complex. This way 
of defining the resultant was introduced by Cayley [GKZ94, app. A], [Wey94], 
[WZ94]. In the particular case in which the complex has just two terms, its 
determinant is nothing but the determinant of the only nonzero differential, 
which is therefore equal to the resultant. In this case, we say that there is a 
determinantal formula for the resultant and the corresponding degree vector m 
is called determinantal. In [WZ94], the multihomogeneous systems for which a 
determinantal formula exists were classified; see also [GKZ94, Sec. 13.2]. Their 
work, though, does not identify completely the corresponding morphisms nor 
the determinantal vectors m, a question we partially undertake. We follow the 
results in [DDOl], which concerned the homogeneous case, inspired also by 
[Jou97]. 

The main result of [SZ94] was to prove that a determinantal formula of 
Sylvester type exists exactly when all defects arc zero. In [SZ94, Thm. 2] 
(recalled in [GKZ94, Sec. 13.2, Prop. 2.2]; see also [McC33, Thm. 4]) aU such 
formulae are characterized by showing a bijection with the permutations of 
{1, ■ ■ ■ ,r} and defined the corresponding degree vector m as in Definition 5.2 
below. This includes all known Sylvester-type formulae, in particular, linear 
systems, systems of 2 univariate polynomials and bihomogeneous systems of 3 
polynomials whose resultant is, respectively, the coefficient determinant, the 
Sylvester resultant and the Dixon resultant. In fact, Sturmfels and Zelevinsky 
characterized all determinantal Cayley-Koszul complexes, which are instances 
of the Weyman complexes when all the higher cohomologies vanish. 

The incremental algorithm for sparse resultant matrices [EC95] relies on the 
determination of a degree vector m. When 5 = 0, it produces optimal Sylvester 
matrices by [SZ94]. For other multihomogeneous systems, [EC95] heuristically 
produces small matrices, yet with no guarantee. For instance, on the system 
of Ex. 5.5, it finds a 1120 x 1200 matrix. The present paper explains the 
behavior of the algorithm, since the latter uses degree vectors following (5.2) 
defined by random permutations. Our results provide immediately the smallest 
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possible matrix. More importantly, the same software constructs all Sylvester- 
type formulae described here. 

Pure Bezout-type formulae were studied in [CKOO] for unmixed systems whose 
support is the direct sum of what they call basis simplices, i.e. the convex 
hull of the origin and another points, each lying on a coordinate axis. This 
includes the case of multihomogeneous systems. They showed that in this case 
a sparse resultant matrix can be constructed from the Bezoutian polynomial, 
to be defined in Section 6, though the corresponding matrix formula is not 
always determinantal. Their Corollary 4.2.1 states that for multihomogeneous 
systems with null defect vector, the Bezout formula becomes determinantal; 
Saxena had proved the special case of all /fc = 1 [Sax97]. In [CKOO, Sec. 4.2] 
they indicate there are r! such formulae and in Section 5 they study bivariate 
systems {n = 2) showing that then, these are the only determinantal formulae. 

Section 6 proves these results in a different manner and characterizes the de- 
terminantal cases for multihomogeneous systems, showing that a null defect 
vector is a sufficient but also necessary condition for a determinantal formula 
of pure Bezout type for any n. Thus, there is an optimal Sylvester-type for- 
mula for the resultant if and only of there is an optimal pure Bezout-type 
formula, (cf. Definition 6.1) This had been proven for arbitrary systems only 
in the bivariate case [CGZ98]. In particular, we explicitly exhibit a choice of 
the differential in the Weyman complex in this case (cf. Theorem 6.13), thus 
partially answering the "challenge to make these maps explicit" of [WZ94, 
Sec. 5.1]. 

Studies exist [CCZ98], [DDOl], [ZhaOO] for dealing with hybrid formulae in- 
cluding Bczout-typc blocks or pure Bezout matrices, and concentrate on the 
computation of such matrices. In particular, [ZhaOO] elaborates on the rela- 
tion of Sylvester and Bezout-type matrices (called Cayley-type there) and the 
transformations that link them. The theoretical setting together with Pfaffian 
formulae for resultants is addressed in [ESOl]. This is made explicit for any 
toric surface in [K02]. In the recent preprint [ACC02], not also the multiho- 
mogeneous resultant but the whole ideal of inertia forms is studied, extending 
results of Jouanolou [Jou80] in the homogeneous case. 

This paper is organized as follows. The next section provides some technical 
facts useful later. Section 3 offers bounds in searching for the smallest possible 
determinantal (hybrid) formulae. Section 4 makes explicit one degree vector 
attached to any determinantal data of dimensions of projective spaces, and dis- 
cusses further techniques for obtaining determinantal formulae. Section 5 and 
6 characterize matrices of pure Sylvester and pure Bezout type respectively. 
In Section 7 we fully describe the formulae of a system of 3 bilinear polyno- 
mials. Then, we provide an explicit example of a hybrid resultant matrix for 
a multidegree for which neither pure Sylvester nor pure Bezout determinantal 
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formulae exist; this example illustrates the possible morphisms that may be 
encountered with multihomogeneous systems. Our Maple implementation is 
described in Section 8. 

A preliminary version of certain results in this paper has appeared in [DE02]. 



2 Preliminary observations 

We consider the r-fold product X := P'^ x . . . x P''" of projective spaces of re- 
spective dimensions li, . . . ,lr over an algebraically closed field of characteristic 
zero, for some natural number r. We denote by n = J2k=i ^fc the dimension of 
X, i.e. the number of affine variables. 

Definition 2.1 Consider d = {di, . . . , d,) G and multihomogeneous poly- 
nomials fo, ■ ■ ■ , fn of degree d. The multihomogeneous resultant is an irre- 
ducible polynomial R{fo, . . . , fn) = R(h,...,ir)Afo, ■ ■ ■ , fn) m the coefficients of 
fo, ■ ■ ■ , fn which vanishes iff the polynomials have a common root in X . 

This is an instance of the sparse resultant [GKZ94]. It may be chosen with 
integer coefficients, and it is uniquely defined up to sign by the requirement 
that it has relatively prime coefficients. The resultant polynomial is itself ho- 
mogeneous in the coefficients of each /j, with degree given by the multihomo- 
geneous Bezout bound (^^ Jdi^ ■■■d'-; [GKZ94, Prop. 13.2.1]. This number 
is also called the m-homogeneous bound [Wam92]. 

Let V be the space of (n -f- 1) tuples / = (/o, . . . , fn) of multihomogeneous 
forms of degree d over X. Given a degree vector m & V there exists a finite 

complex K. = K.{m) of free modules over the ring of polynomial functions 
on V [WZ94], whose terms depend only on {li, . . . ,lr),d and m and whose 
differentials are polynomials on V satisfying: 

(i) For every given / we can specialize the differentials in K. by evaluating at 
/ to get a complex of finite-dimensional vector spaces, 
(n) This complex is exact iff R{fo, . . . , /n) 7^ 0. 

In order to describe the terms in these complexes some facts from cohomology 
theory are necessary; see [Har77] for details. Given a degree vector m e Z'', 
define, for 1/ e {—n, . . . ,n-\- 1}, 

Ku{m)= HP-''{X,m-pd)("p'\ (1) 

pe{0,...,n+l} 

where for an integer r-tuple m', H'^{X, m') denotes the q-th cohomology of X 
with coefficients in the sheaf 0{m') such that its global sections H^{X,m') 
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are identified with multihomogeneous polynomials of (multi)degree m' . By the 
Kiinneth formula, we have 

jfcG{0,/fc} r 

H'i{X, m-pd)^ (g) W'iW'^mu - pd^), 

iiH \-jr=q k=l 

where q = p — u and the second sum runs over all integer sums Ji + ■ ■ ■ + 
jr = Q,jk £ {0,/fe}. In particular, {¥'■'', a^) is the space of all homogeneous 
polynomials in + 1 variables with total degree 0;^. By Serre's duality, for 
any a e Z^, we also know that 

H'^iX, a) ~ H'^-^X, {-h - 1, 1) - ay, (2) 

where * denotes dual. We recall Bott's formulae for these cohomologies. 

Proposition 2.2 For any m G Z*", H^>'(P^>' ,mk—pdk) = rrik—pdk > —Ik, 
H°{¥^'',mk - pdk) = <^ mfe - pdk < 0, for k E {1, . . . , r}. Moreover, 



■mk + pdk - 1\ 
Ik I 
rrik - pdk + lk\ 



dimii■''=(P'^mfe -p4) = 
dimH'^{¥\mk-pdk) ^ 



Consequently, 

ifce{o,«fc} r 

(limH'i{X,m-pd) = ^ [] dimi^^'=(P^^ - p4), 

JiH 1->=9 fe=l 

and 

dim ii'^(m)= ^ r ^ M dimii'^'-'^(X,m-p(i). 

pG[0,n+l] \ P J 

Definition 2.3 Given r and (/i, . . . , Z^), (rfi, . . . , 4) G , define the defect 
vector 6 e Z'' (just as m [SZ94], [WZ94]) by 6k := h - f^l- Clearly, this 
is a non-negative vector. We also define the critical degree vector p & W by 
Pk := {n + 1)4 - Ik - I, for all k ^ 1, . . .r. 

Lemma 2.4 [WZ94] For any i e [r] := {1, . . . , r}, dik < di + k <^ 6i = ^ 
mm{li,di} = 1. 

Let us establish a general technical lemma. 

Lemma 2.5 For any /c e {1, . . . , r}, < (Z^ — 5fc)4 — ^fc < 4 — 1- 
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Proof. By definition, 5k = h ~ \h/dk] ^ \lk/dk] ^ h - h ^ {lk + tk)/dk for 
some integer tk sucli tliat < tk < dk — 1- □ 

We detail now tlie main results in [WZ94]. They show (Lemma 3.3(a)) that 
a vector m G Z*" is determinantal iff K_i{m) = K2{m) = 0. They also prove 
in Theorem 3.1 that a determinantal vector m exists iS Sk < 2 for all A; e [r]. 
To describe a differential in the complex from K^{'m) to K^+i{m), one needs 
to describe all the morphisms 5py from the summand corresponding to an 
integer p to the summand corresponding to another integer p', where both 
p,p' E {0, . . . ,n + 1}. [WZ94, Prop. 2.5, 2.6] proves this map is when p < p' 
and that, roughly speaking, it corresponds to a Sylvester map {go, . . . , gn) — > 
J27=o 9ifi when p — p' thus having all nonzero entries in the corresponding 
matrix given by coefficients of /o, . . . , /«. For p > p' + 1, the maps are 
called higher-order differentials. By degree reasons, they cannot be given by 
Sylvester matrices. Theorem 2.10 also gives an explicit theoretical construction 
of the higher-order differentials in the pure Bezout case (cf. Definition 6.1). 



3 Bounds for determinantal degree vectors 



This section addresses the computational problem of enumerating all deter- 
minantal degree vectors m Ell . The "procedure" of [WZ94, Sec. 3] "is quite 
explicit but it seems that there is no nice way to parametrize these vectors" , as 
stated in that paper. Instead, we bound the range of m to implement a com- 
puter search for them. In Section 4 we will give an explicit choice of degree 
vector m for each determinantal data ( l\, . . . ,lr\d\, . . . ,dr\ 

Given k e {1, . . . , r} and a vector m eU define as in [WZ94]: 
"fe \^) = {p E Z : — < p < 



dk d_ 



k 



Let Pk{m) be the real interval (^,^^^^], so Pk{m) = Pk{m) n Z. Using 
Lemma 3.3 in [WZ94], it is easy to give bounds for all determinantal vectors 
m for which all Pkim) ^ 0. 

Lemma 3.1 For a determinantal m G Z*^' and for all k E {1, . . . ,r}, Pk{iTi) 7^ 
implies, max{— 0?^, —Ik} < ruk < dk{n -|- 1) — 1 -|- m.m{dk — Ik, 0}. 

Proof. By [WZ94, Lem. 3.3(b)], Pfc(m) C [0,n+ 1] ruk/dk > -1 and 
{mk + lk)/dk > 0, which imply the lower bound. Also, {mk + lk)/dk <n + 2 ^ 
rrik < {n + 2)dk — Ik — I and rrik/dk < n + 1 rrik < {n + l)dk — 1 yield the 
upper bound. Notice that the possible values for nik form a non-empty set, 
since the two bounds are negative and positive respectively. □ 
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Now, p e Pfc(m) iff H^^iF^^ruk - pdk) = H'^{F^'',mk - pdk) = 0. Thus, a first 
guess could be that all determinantal vectors give Pk{in) ^ 0. But this is not 
the the following example shows: 

Example 3.2 Set I — (l,2),d = (2,3). We focus on degree vectors of the 
form m = (2/ii,3/i2), for /ii,/i2 € Z. Then, for all such m, the sets Pk{m), for 
k = 1,2, are empty. Nevertheless, there exist 4 determinantal vectors of this 
form, namely m = (4, 3), (0, 6), (2, 6) or (6, 3). Moreover, the vector m = (6, 3) 
gives a determinantal formula with a matrix of size 88, which is closer to the 
smallest possible one which has size 72. The largest determinantal formula is 
given by the determinant of a square matrix of size 180. Note that the degree 
of the multihomogeneous resultant is 216. More details on this example are 
provided in Section 8. □ 

We wish now to get a bound for those determinantal vectors for which some 
Pfc(m) is empty. Let [-jk G {0, 1, . . . , — 1} denote remainder after division 
by 4. 

Definition 3.3 Given m e and /c e {1, . . . , r}, define new vectors m' , m" e 
7/ whose j-th coordinates equal those of m for all j ^ k and such that 
m'k ^ruk + dk- [mk]k - 1 > ruk > m'l ^ ruk - [m^ + lk\k- 

Lemma 3.4 The vectors m',m" differ from m at their k-th coordinate if 
Pk{m) = 0. 

Proof. Let us write = jdk + {mk]k for j e Z. Then, ^ = j + I^^. 
If Pfc(m) = 0, < J + 1 ^ jdk + [mk]k + lk < {j + 1)4, so [mk]k < 

dk — h — ^ =^ h ^ — [^fc]fc + 4 ~ 1 and thus 1 < 4 ~ [^^fc]fe ~ 1- Also, 
[rrik + lk]k > 1 because [ruk + lk]k = ^ (m^ + 4)/4 e Pk{m). □ 

Lemma 3.5 If m E U' with Pk{rn) = and H'^ {¥'■>', nik — pdk) = (resp. 
/^"''=(P'^mfe -p4) = O;, then Pk{m') 7^ and H^{F^'^,m'^ - pdk) = (resp. 
H'''' {¥'''', m'f.— pdk) — 0), where m'j^ — THk + dk — [mk]k — i as in Definition 3.3. 

Proof. Write nik = jdk + [mk]k for some integer j G Z. To prove Pk{m') ^ 
we show j + 1 e Pk{m'), i.e. ^ < j + 1 < ^^^^ < (j + 1)4 <m'f^ + lk 

^ jdk + 4 ~ 1 < (j + ^)dk < jdk + dk — 1 + Ik-i which is clearly true since 
1 <lk- Now, if''=(P'*=, ruk — pdk) = <^ + > p4 hence m'^ + /fc > p4 
because m'^ > rUk- //'^(P''-, — pdk) = -v^ < pdk so j < p — 1. Hence, 
rUk — [rriklk = jdk < (p — 1)4 f^k — [f^k\k + dk < pdk which is the desired 
conclusion. By [WZ94], Pk{m') C [0, n + 1] from which j G {-1, 0, . . . , /}. □ 

Lemma 3.6 If m G Z** with Pkijn) = and H^{¥'-'',mk — pdk) = (resp. 
H^k(^ik^mk- pdk) = 0), then Pk{m") 7^ and H^{¥^\ml - pdk) = (resp. 
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Hh(f'ik^ — pdk) — 0), where m'l — ruk — [ruk + 4]fe o-s in Definition 3.3. 

Proof. Write rrik + h = jdk + [rrik + lk]k for some integer j > 0. To prove 
Pk{m';:) ^ we show it contains j, i.e., ""^-^2^^"^'' < 3 < ""'"'""'d;^''''^'' ^ 
"rrik - [rUk + lk]k < jdk < mk - [rrik + lk]k + k jdk - k < jdk < jdk, which 
is clearly true since < Ik- H^iF'-'^ ,mk — pdk) = -v^ rrik < pdk ^ fn'l < pdk 
because m'/, < rrik, hence //"(P'*^, m'^' — pdk) = 0. i/''=(P''=, — pdk) = <;=^ 
f^k + h ^ pdk, then j > p. Hence rrik + h — [m-k + h]k > pdk which finishes 
the proof. By [WZ94], Pfc(m") C [0, n + 1] hence j G {0, . . . , n + 1}. □ 

Lemmas 3.5 and 3.6 imply 

Theorem 3.7 For any determinantal m G Z*", define vectors m',m" G 1/ as 
in Definition 3.3 which differ from m only at the k-th coordinates, 1 < k < r, 
such that Pk{m) = 0. Then Pk{m') ^ i,Pk{m") ^ and hath m',m" are 
determinantal. 

Corollary 3.8 For a determinantal m G Z*^' with Pk{fn) = for some k G 
{1, . . . , r}, we have < rrik < dk{n + 1) — 4 — 1. 

Proof. Since m^, rri'l define Pkim') ^ 0, Pk{m") ^ 0, we can apply Lemma 3.1. 
We use the lower bound with m'^ because m'^ < rrik < '^'k- Pkijn) = ^ => dk > 
Ik, so ml = mk- [rrik + lk]k> -h ^ rrik > [rrik + lk]k - h > 'i- - k, because 
[rrik + ^ 1 by the proof of Lemma 3.4. If rrik < 0, for Pk{m) to be empty 
we need ruk + /fc < ruk < —Ik which contradicts the derived lower bound; 
so ruk > 0. For the upper bound, = rrik + dk — [mk]k — 1 < dk{n + 1) — 1 =^ 
mk < dkU + [mk]k < dk{n + 1) ^ Ik — 1; the latter follows from [mk]k < dk — lk, 
(Lemma 3.4). {rrik + lk) / dk <n+l — {l/dk) < n+ 1 implies the inclusion of the 
half-open interval in (0, n + 1). The possible values for form a non-empty 
set, since the lower bound is zero and the upper bound is dk{n -\- 1) — Ik — 1 > 
dk — I > since dk > h > D So, in fact, the real interval Pk C (0, n -|- 1). 

Corollary 3.9 For a determinantal m G and k G [r], max{— d^, — Z^} < 
mk < dk{n + 1) - 1 + mm{dk - h, 0}. 

This implies there is a finite number of vectors to be tested in order to enu- 
merate all possible determinantal m. This could also be deduced from the fact 
that the dimension of KQ^m) equals the degree of the resultant. Corollary 3.9 
gives a precise bound for the box in which to search algorithmically for all de- 
terminantal m, including those that are "pure" in the terminology of [WZ94]. 
Our Maple implementation, along with examples, is presented in Section 8. 

If we take r = 1 = /i and m = 2di — 1 we obtain the classical Sylvester formula 
and the upper bound given by Corollary 3.8 is attained. At the lower bound 
m = -1, the corresponding complex H^{¥^, -1 - 2d) = H°{¥^,2di - 1)* 
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H^(P^,di — ly = {H^{¥^,di — 1)*)^ yields the same matrix transposed. In 
addition, the system of 3 bilinear polynomials in Section 7.1 admits a pure 
Sylvester formulae with m = (2,-1), which attains both lower and upper 
bounds of the Corollary. The bounds in Lemma 3.1 can also be attained (see 
Ex. 5.5 continued in Section 8) hence Corollary 3.9 is tight. It is possible that 
some combination of the coordinates of m restricts the search space. 



4 Explicit determinantal degree vectors 



We focus on the case of 5^ < 2 for all k, which is a necessary and sufficient con- 
dition for the existence of a determinantal complex [WZ94] . We shall describe 
specific degree vectors m that yield determinantal complexes. 

Theorem 4.1 Suppose that < 5^ < 2 for all k = 1, . . . , r and tt : [r] ^ [r] 
is any permutation. Then, the degree vector rn" e 1/ with 

ml^ \ l-5k+ E h\dk- k, (3) 

\ ir{j)>n{k) J 

for k — 1, . . . ,r defines a determinantal complex. 

Proof. First, Ki ^ and Ko ^ because they each contain at least one 
nonzero direct summand, namely H^{X,m^ — pd) for p = 1, respectively. 
To see this, it suffices to prove m% — dk > 0, k — 1, . . . ,r, which follows from 
{1 — 6k + lk)dk — h — dk > <^ {—Sk + lk)dk — ^a; > 0. This holds by Lemma 
2.5. 

To demonstrate that i^^2 = we shall see that every one of its direct summands 
H'^{X, rri^ — {q + 2)d) vanishes for q — J2-K{j)ej hi where J is any proper subset 
of [r]. Recall that the case J — [r] is irrelevant by the definitions of Section 2. 
It is enough to show i/°(P''=,m^ ~ + '^)dk) = 0, for some k with 7i{k) ^ J. 
Let k be the maximum of the indices verifying 7r{k) ^ J; such a k G [r] always 
exists because J 7^ [r]. Since (1 — + h)dk — h< '^dk by Lemma 2.5, it easily 
follows that 

ml<{2+ h)dk. (4) 

7r(j)eJ 

It now suffices to establish K_i{m) — 0. Since this module has no zero co- 
homology summand, let q = Y^Tr{j)&jhi with J ^ ^. Let k & J such that 
7r(j) > iT{k) for all j G J. Then, 6k < 2 implies (2 — 6k)dk > 0, from which 
^k~il~^)dk > — /fe, and so H'^{X,m'' — {q—l)d) = for any direct summand 
ofK_i{m). □ 
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When the defects are at most 1, we can give another exphcit choice of deter- 
minantal degree vector for each permutation of [r]. 

Theorem 4.2 Suppose that < 6k < 1 for all k — 1, . . . ,r and tt : [r] — > [r] 
is any permutation. Then, the degree vector rri^ e Z'' with 



\ 7r(i)>7r(A;) / 

for k — 1, . . . ,r defines a determinantal complex. 

Proof. Let us modify the previous proof. First, Kq because it contains 
H^{X,rn^ -pd) for p = 0. This follows from {-Sk + lk)dk - k > ^ 
{—^k + h)dk — ^fc > for all k. This holds by Lemma 2.5. 

Next, Ki because for p = n + 1, if'*(P''=,m^ — pdk) 7^ for all k. This 
follows from {—5k + n)dk — Ik < {n + l)dk — Ik —Sk < 1. 

To demonstrate that K2 = we repeat the corresponding argument in the 
proof of Theorem 4.L Then, it suffices to show that {—Sk + lk)dk — h < '^dk, 
which is weaker than the inequahty before (4) above. Similarly, to estabhsh 
K_i{m) — the steps in the previous proof lead us to showing m^ — (g— > 
—Ik <= {1 — Sk)dk > which holds for 5^ < 1. □ 

One can easily verify that the vectors in both Theorems above satisfy the 
bounds of Corollary 3.9. 

A natural question is how to give an explicit degree vector yielding a deter- 
minantal formula of smallest size. When r — 1, this is the content of Lemma 
5.3 in [DDOl]. Even for r — 2, this seems to be a difficult task in general. 

Example 4.3 Set Z = (2, 2), d = (3, 2). Let tti : [2] [2] denote the identity 
and 772 '■ [2] — >■ [2] the permutation which interchanges 1 and 2. The two 
determinantal vectors defined according to (3) are m^^ = (10,2) and 777,7^2 = 
(4, 6), which yield determinantal matrices of respective sizes 396 and 420. The 
two determinantal vectors defined according to (5) are m^j = (7, 0) and mTra — 
(1, 4), which yield determinantal matrices of respective sizes 756 and 780. But 
as we will sec in §6, these latter formulae give instead the vectors providing 
the smallest determinantal formulas when all defects are 0. These vectors 
can be computed by the function comp_m of our implementation, discussed in 
Section 8. 

Also, when r = 1 it is shown in [DDOl] that the smallest formula is attained 
for "central" determinantal degree vectors. In this example, all degree vectors 
m satisfy that either 1 < ?77i < 5 and 4<m2 <7or7<mi < 11 and 
< 77i2 < 4; these bounds are computed by the routines described in §8. The 
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smallest resultant matrix has size 340 x 340, and corresponds to the degree 
vectors (3, 6) or (9, 1), which are in a "central" position among determinantal 
degree vectors, in other words, their coordinates lie in a "central" position 
between the respective coordinates of other determinantal vectors. However, 
the vectors {k, 4) are determinantal for k from 1 to 5, but the size of a matrix 
corresponding to the vector (3, 4) (both of whose coordinates lie between the 
coordinates of the vectors (1,4) and (5, 4)), equals 580. So, it is bigger than the 
size of a resultant matrix associated to the vector (5, 4), which is of dimension 
540. □ 



In the above example, the two degree vectors giving the smallest resultant 
matrices satisfy (9, 1) + (3, 6) = (12, 7), which is the critical vector from Defi- 
nition 2.3. It is clear that for any determinantal m, the vector (12, 7) — m is 
also determinantal yielding the same matrix dimension. This is a consequence 
of Serre's duality recalled in (2). The general statement in summarized in the 
next Proposition. 

Proposition 4.4 Assume m, m' e satisfy m + m' = p, the latter being the 
critical degree vector of Definition 2.3. Then, Ky{m) is dual to Ki_y{m') for 
all V In particular, m is determinantal if and only ifm' is determinantal, 
yielding matrices of the same size, namely <\mi{KQ{m)) — <\mi{Ki{m')) . 

Proof. Based on the equality m + m' = p we deduce that for all p = 0, . . . , n + 
1, it holds that (m' — pd) — {—h — 1, . . . , —Ir — 1) — (m — (n + 1 — p)d). 
Therefore, for all g = 0, . . . , n, Serre's duality (2) implies that H'^{X, m! — pd) 
and //"""^(X, m — (n+ 1 —p)d) are dual. Since {n + 1 — p) — {n — q) = l — {p — q), 
we deduce that K^{m) is dual to A'i_i^(m') for all z/ G Z, as desired. Observe 
that in particular i^_i(m) ~ K2{m')* and Ko^m) ~ Ki{m')*, the latter giving 
the matrix dimension in the case of determinantal formulae. □ 



We end this section making explicit a consequence of Proposition 3.7 in [WZ94] 

(and giving an independent proof), which gives a generalization of the char- 
acterization of determinantal complexes in [SZ94], [WZ94]. This result will 
allow us to give in §6 explicit expressions for all degree vectors yielding a 
determinantal formula of smallest size when all defects vanish. 

Theorem 4.5 There exists a determinantal vector m such that the Wey- 
man complex is reduced to only one non-zero cohomology group on each of 
Ko^m), Ki{m) if and only if all defects vanish, i.e. Sk = for all k = 1, . . . ,r. 

Proof. Recall that for each p G {0, . . . ,n-\-l} there exists at most one integer 
j such that {X, m — pd) 7^ in Ky{m) where p = j — v [WZ94, Prop. 2.4]. 
In fact, let A{p) := {k : — pdk < —h} and B{p) := {k : rrik — pd^ > 0}. 
Denote j{p) :— J2keAip)h- Then H^' {X,m — pd) — for all j' ^ j{p) and 
H^^^ {X, m - pd) ^ iff A{p) U B{p) = {1, . . . , r}. 
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The assumption of the theorem means that there exist exactly two integers 
Pi,P2 e {0, . . . , n + 1} for which H^'^^'^ {X,m — pid) ^ 0,i = 1,2; cf. also 
[WZ94, Lem. 3.3(a)]. Then, for any p G {0, . . . ,n + 1} \ {pi,P2}, there exists 
k such that k ^ Ai^p) U B{p), i.e., p G Pfc(m). The latter set is defined at the 
beginning of Section 3. Then, {0, . . . , n + 1} \ {pi,P2} ^ U^^iPfc(m) and so 

l^ + ,,, + l^^ri<j2i^Pk(m)<j2\^], 

k=l k=l "'k 

where the first inequality uses the fact that # U^^^ -Pfe(^) ^ Z)fe=i#-Pfc(^) 
and the second follows from the definition of Pk{m). Since Ik > [^] for all k, 
we deduce that Ik = \^~\ , and this can only happen iff Zjt = 1 or dk — 1- □ 



5 Pure Sylvester-type formulae 



This section constructs rectangular matrices of pure Sylvester-type that have 
at least one maximal minor which is a nontrivial multiple of the sparse resul- 
tant, coming from a complex of the form: 

> K2{m) Ki{m) Ko{m) K_i{m) = 0, 

where Ki{m) = H^{X, m — dY~^^, and KQ{m) — H^{X, m) for a non negative 
vector m G Z>q. 

We assume that H^{X,m — pd) 7^ for z/ = 0, 1 and p — = 0. This im- 
plies i7°(P'*,mfc — udk) 7^ for all k G {1, . . . ,r}. Moreover, we must have 
HP'^[X, m — pd) — ior p — u — Z^jg j Ij where J is any subset satisfying 
07^Jc{l,...,r}. Note that we do not require in general that K2{m) ^ 0. 

Lemma 5.1 If m ^ m' & Z>q yield a Sylvester-type matrix and m'f. > ruk for 
a// A; G {1, . . . , r} then, the Sylvester matrix associated to m' is strictly larger 
than the Sylvester matrix associated to m. 

Proof. We must show dimii'o(m') > dimi^o(^)) i-G- dimif°(P''=, m'^ — pd^) > 
dimH\¥\mk-pdk) for A; G {1, . . . , r}, i.e. {^'>.-fj.+h\^ > ^u-pd,+i,\^^ 
cohomology is nonzero, thus m'^, — pd^ > nik — pdk > and this implies the 
desired inequality because = {s + Ik) ■ ■ ■ {s + l)//fc!- The inequality is 

strict since there exists an index k such that m'^ > ruk- □ 

Definition 5.2 For each choice of a permutation tt : {1, . . . , r} ^ {1, . . . , r}, 
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consider the degree vector rri^ defined by 

"^fc-= h-+ Yl lj\dk-lk, k^l,...,r. 

\ 7r(i)>7r(fc) / 

When all defects are zero, these are the vectors defined in [SZ94] yielding 
determinant al Sylvester formulae and they also coincide with those defined in 
(3) in Theorem 4.1. 

Lemma 5.3 If m E Z>q yields a Sylvester-type matrix, it is possible to define 
a permutation tt : [r] — > [r] such that for i = 7r^^(l) it holds that rui > mj . 
Moreover H^{F''\mi —pdi) ^ where p < 1 + I]7r(j)Gj ^i? for any subset J such 
that^^ JC {2, ...,r}. 

Proof. For p = n+l,z/=l, a necessary condition is that H"'{X,m — (n + 
l)d) = 0. Hence, there exists i G [r] : H'-'(¥'-\mi — (n + l)di) = mj — (n + 
i)di > —li rUi > by choosing 7r(i) — 1. For any p as in the statement, 
if°(P'%mi — pdi) 7^ <^ mi > pdi. Since rrii — {n + l)di > —k it suffices to 
prove (n + l)di — U > di{l + Y^j^ilj) > dip. The latter inequality is obvious 
for all p, whereas the former reduces to /jdj > li which holds since di> 1. □ 



Theorem 5.4 A degree vector m G Z>g gives a Sylvester-type matrix iff there 
exists a permutation tt such that rrij > mj for j = 1, . . . ,r. Moreover, the 
smallest Sylvester matrix is attained among the vectors rn^ . 

Proof. We prove the forward direction by induction on k = l,...,r. As- 
sume m gives a Sylvester-type complex and consider the necessary condition 
Ki(m) — H^{X, m — d). The base case k — 1 was proven in Lemma 5.3. The 
inductive hypothesis for k G {l,...,r — 1} specifies which cohomologies vanish 
and which not, where m„ > m^, tt{u) < k. In particular, for all subsets J such 
that 7^ J C {1, . . . , r} \ {1, . . . , /c}, p=l + Y.-!T(j}eJ Po=P + ^v, for some 
V such that n{v) < k, we assume: 

H'- (P^" , m„ - pod^) = 0, (P'" , m„ - pd^) ^0. (6) 

For the inductive step, we exploit the necessary condition that H^~^[X,m — 
pd) = forp = 1+Zl7r(j)>fc By (6), there exists i such that (^P'% — pd^ 
0. Then, rUi > pdi—li — where we define n{i) = k-\-l. To complete the step, 
we show (¥'■", ruu — pd^ ^ where 'k{u) < k -\- 1, p — 1 -\- Ij^ and 

any subset J such that 7^ J C {A;-|-2, . . . , r}. The non- vanishing of the coho- 
mology is equivalent to m„ > pdu- It suffices to prove mj > di{l-\-J2Tr{j)>k+i h)- 
By definition, this reduces to — + diU > -v^ lij > 1. The converse direction 
follows from analogous arguments as above. The claim on minimality follows 
from Lemma 5.1. □ 
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This gives an algorithm for finding the minimal Sylvester formulae by testing 
at most r! vectors m'^ , which is implemented in MAPLE (Section 8). To actually 
obtain the square submatrix whose determinant is divisible by the sparse 
resultant, it suffices to execute a rank test. These matrices exhibit quasi- 
Toeplitz structure, implying that asymptotic complexity is quasi-quadratic in 
the matrix dimension [EP02]. Observe that Pk{m'") 7^ because there exists 
p e Z such that p — T,'K(j)<n{k) h such that ml < dkP — ml + Ik for all k. 

Example 5.5 Let I — (2,1,1), d = (2, 2, 2); the degree of the resultant is 960. 
Let (7 = 7r~^ be the permutation inverse to tt; then the corresponding degree 
vector can be written as m^^^^ := (l + J2j>k^cr(j)j d„{k) — la(k)- Here is a list 
of the 6 = 3! degree vectors m^, among which we find the smallest Sylvester 
matrix of row dimension 1080, whereas the sparse resultant's degree is 960. 
Also shown are the permutations a and the corresponding matrix dimensions. 
The symmetry between the last two polynomials makes certain dimensions 
appear twice. 



(8,5,3) 


a = (1,2,3) 


1080 X 1120 


(8,3,5) 


(1,3,2) 


1080 X 1120 


(6,9,3) 


(2,1,3) 


1120 X 1200 


(4,9,7) 


(2,3,1) 


1200 X 1440 


(6,3,9) 


(3,1,2) 


1120 X 1200 


(4,7,9) 


(3,2,1) 


1200 X 1440 



Our Maple program, discussed in Section 8, enumerates 81 purely rectangular 
Sylvester matrices (none of which is determinantal) . All Sylvester matrices not 
shown here have dimensions 1260 x 1400 or larger. □ 

The map Ki{m) KQ^ni) is surjective, i.e., the matrix has at least as 
many columns as rows. In searching for a minimal formula, we should reduce 
dimKo{m), i.e., the number of rows, since this defines the degree of the ex- 
traneous factor in the determinant. It is an open question whether dimKo{m) 
reduces iff diinKi{m) reduces. In certain system solving applications, the ex- 
traneous factor simply leads to a superset of the common isolated roots, so 
it poses no limitation. Even if it vanishes identically, perturbation techniques 
yield a nontrivial projection operator [DEOl]. 

It is possible to obtain a pure Sylvester-matrix whose determinant equals the 
multihomogeneous resultant when the complex has as only non-zero terms 
Ki{m) = W{X,m- {j + l)d)^^+'\ and Ko{m) = W{X,m - jd)^"^'^ for any 
j — 0, . . . , n. But in this case we deduce from Theorem 4.5 that all defects 
vanish. So (cf. [SZ94]) there exists a pure Sylvester-type determinantal formula 
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associated to a non negative degree vector (i.e. nik > 0, k = l,...,r), or 
equivalently. for j = 0. Thus, if a determinantal data ( /i, . . . , Z^; di, . . . , dj.) 
admits a degree vector yielding a two term Sylvester complex for some j, it 
admits such a formula for j — as well. Hence, concentrating on j = is not 
restrictive in the case of determinantal complexes. The Sylvester-type matrices 
for positive j correspond to degree vectors with some negative entries, unlike 
the assumption in [SZ94, p. 118]. We show such an example in the bihnear case 
in §7. Thus, the first part of conjecture [SZ94, Conj.3] can be true only for 
non-negative degree vectors. 



6 Pure Bezout-type formulae 

In this section, we will study the following complexes: 

Definition 6.1 A Weyman complex is of pure Bezout type if K-i[m) = 0, 
Ki{m) = H^^+-+^r (x,m - {n + l)d) and Ko{m) = H^m). 

Weyman complexes of pure Bezout type correspond to generically surjective 
maps 

jjh+-+ir m-{n+ l)d) H\m) (7) 

such that any maximal minor is a nontrivial multiple of the multihomogeneous 
resultant. In fact, we shall show that the only possible such formulae are deter- 
minantal (i.e. K2{m) — 0). We shall exhibit the corresponding differential in 
terms of the Bezoutian and characterize the possible degree vectors. We show 
that there exists a pure Bezout-type formula iff there exists a pure Sylvester 
formula. We remark that the dimension of the matrix with pure Bezout coef- 
ficients equals the dimension of the Sylvester matrix divided by n -|- 1. Now 
we can generalize results in [CKOO], [Sax97] (cf. Section 1). 

Theorem 6.2 There exists a determinantal formula of pure Bezout type iff 
for all k either 1^ = 1 or d^ = 1, i.e. all defects vanish. 

Proof. This is just a special instance of Theorem 4.5, where the only non-zero 
cohomologies in the complex correspond io pi — Q,p2 — n -\- 1. □ 

Let us study degree vectors yielding pure Bezout formulae, which will then 
provide the smallest determinantal formulae in case all defects vanish. 

Definition 6.3 For each choice of a permutation tt : {1, . . . , r} ^ {1, . . . , r}, 



16 



let us define a degree vector 



iK ■■= -Ik + 4 ij, k = l,...,r. 

T(j)>7r(fe) 

When all 6k = these are precisely the vectors defined in (5) in Theorem 
4.2. Note that our assumptions in (7) imply H^^ (F^^rrij — (n + l)dj^ ^ 0, 

(F^^TJlj) 0. 

Lemma 6.4 The existence of any pure Bezout formula implies < rrij < 
{n + l)dj — Ij, for all j. 

In fact, the of Def. 6.3 satisfy these constraints for all permutations tt. 

Proof. The non negativity of all rrij is deduced from the fact that H'^{X, m) ^ 
0. On the other side, the non vanishing of H"'[X,m — {n + l)d) implies the 
other inequality. □ 

Lemma 6.5 IfmEZ^ yields a pure Bezout-type complex, then there exists a 
permutation tt : [r] — > [r] such that i — '^'^{1) verifies rrii > rn^ and 

H^' (P^Smi -{q + li + y)d,) =0, z/ = 0, -1, 
for any J C{1,.. .,r}\{i}, J 

Proof. Since H"'{X,m — nd) = 0, there exists an index i G {1, . . . ,r} such 
that rrii — ndi > — /j. It is enough to define 7r(i) = 1. □ 

Theorem 6.6 If m E yields a pure Bezout-type complex, it is possible to 
find a permutation tt such that the degree vector m verifies rrii ^ /'^'^ 
i = l,...,r. 

Proof. We use induction; the base case follows from Lemma 6.5. The in- 
ductive hypothesis, for /c e {!,..., r — 1}, is: there exists a subset U C 
{!,..., r}, \U\ — k, such that Tr{u) < k, m„ > for all w e C/ and 

H'- (P'", m„ -{q + lu + y)du) =0, z/ = 0, -1, (8) 
i/°(P^",m„-gd„) 7^0, q = Y.hi 
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for all J C {1, . . . , r}\f/, J 7^ 0. Now the inductive step: The hypothesis on Kq 
implies Hp{X, m—pd) = for p = J2j^u h- Considering the inequality in (8) for 
q = Pi there exists i G [r] \ f/ such that H^^ — pdij = ^ mi + li > pdi 

i.e. rrii > mj for 7r(i) = k + 1 because j ^ U TT{j) > k + 1. It suffices 
now to extend (8) for q' = EjeJ'h where (/} ^ J' G [r] \ {U U {«}). First, 
TTii+li > pdi > {q'+li)di implies the equations below. Second, > —k+pdi — 
{p — lijdi -\- li{di — 1) > q'di yields the inequality, so 



H'^ (pSm^ - {q' + k + u)d,) =0, u = 0,-l. (9) 

(F'\mi-q'di) j^O. 

Then, m satisfies the hypothesis K_i{m) = for z/ = — 1 because every 
summand in K_i{m) contains some cohomology as in (9). Since p > ^ q — 
p — v > 1 no summand has only zero cohomologies. By Lemma 6.4 and (9) 
for u — 0, m gives Ko{m) — H^{X, m) because H^^{¥''%mi — pdi) — for any 
p,i. □ 



Theorem 6.7 A pure Bezout and generically surjective formula exists for 
some vector m iff it equals rn^ of Definition 6.3, for some permutation tt, and 
all defects are zero. 

Proof. It suffices to consider Ki{m) = H"'(X, m — {n + it is nonzero by 
Lemma 6.4. We prove by induction that m = m'^ by using the fact that all 
other summands in (1) for Ki{m) vanish. For H'^[X,m — d) to vanish, there 
must exist i e [r] such that iJ°(P'% mj — di) = Q ^ rrii < di. Hence we need to 
define 7r(i) = r because 7r(i) < r =^ ml > —li + di{li + l) — di + li{di — l) > di. 
Moreover, mJ — —li + dik < di ^ Si — hy Lemma 2.4. 

There is a unique integer in [m^, di) because + 1 > di 4^ —li + dili + 1 > 
di {h — l){di — 1) > 0. Hence rrii — di — 1 < di{q + 1) for any q > 0, 
therefore i?°(P'%mi — (g + l)di) — 0. Furthermore, for q > li, if'*(P'%mj — 
(1 + q)di) rrii + li < {q + ^)di k — 1 < qdi which holds. This proves 
the inductive basis. The inductive hypothesis is: for all u & U C [r], where 
1 6'"! = k, 7t{u) > r — k, then Su = 0, m„ = and 

i7°(P^",m„ - (1 + Yl h)du) = ^ i/'"(P'",m„ - (1 + Y.h)du)i (10) 

7r(j)<7r(w) j&J 

for all J such that U <Z J <Z [r\. For the inductive step, consider that 
H'^{X,m — (1 + q)d) must vanish for q — Y^j^uh- None of its summand 
cohomologies i?'" (^P'",m„ — (1 + q)dy^ vanish due to the last inequality. So 
there exists i such that if°(P'% — (1 + q)di) = <^ < (1 + q)di. 



18 



Hence -n{i) = r — k so that rrii = = —k + diJ27j:{j)>r-kh < (1 + 
—li + dili < di Si = by Lemma 2.4. No larger mj works because is 
the maximum integer strictly smaller than (1 + q)di. And 7r(i) < r — k would 
make too large. Now extend the inequaUty (10) to J' where {U U {i}) C J' 
and observe < di J2-,r{j)>r-k h < di{l + Ejgj' ^j)- 

The hypothesis is proven for all U C [r], including the case \U\ = r. For 
the converse, assume there exists a permutation tt such that m = and 
all defects vanish. Then KQ{m), Ki{m) satisfy all conditions for a pure Bezout 
formula. Furthermore, ir_i(m) = 0, hence the formula is generically surjective. 
□ 

The condition K2{m) = 0, which yields a square matrix, is obtained by the 
hypothesis of a pure Bezout and generically surjective formula; i.e., there is 
no rectangular surjective pure Bezout formula. 

Corollary 6.8 If a generically surjective formula is of pure Bezout type, then 
it is determinantal. Furthermore, for any permutation tt, the matrix is of the 
same dimension, i.e. dimi^o(^) = degi?/(n + 1). 

6.1 Explicit Bezout-type formulae 

We start by defining a dual permutation tt' to any permutation tt. 

Definition 6.9 For any permutation tt : [r] ^ [r], define a new permutation 
tt' : [r] — >• [r] by 7r'(i) = r + 1 — 7r(i). 

Lemma 6.10 Assuming all defects are zero, + m"^' = p for any permuta- 
tion TT : [r] — > [r], where p &W is the critical vector of Definition 2.3. 

Proof, + = di{n + /,) — 2/j for all i because the sum in the parenthesis 
includes {Ij : 7r(j) > 7r(i)} U {Ij : > 7r'(i)}, and latter set is {Ij : 7r{j) < 
7r(i)}. So + mf = dik - k + pi - d^ + 1 = pi + di{li - 1) - {k - 1) = 
Pi + {k — ^){di — I) — Pi because of the zero defects. □ 

Denote by Xi (resp. Xij) the i-th variable group (resp. the j-th variable in the 
group), i e [r],j = 0, . . . , Z^. Introduce r new groups of variables yi with the 
same cardinalities and denote by yij their variables. 

Given a permutation tt, let the associated Bezoutian be the polynomial B'"{x, y) 
obtained as follows: First dehomogenize the polynomials by setting Xiq = 
1, i = 1, . . . , r; the obtained polynomials are denoted by /o, . . . , fn- Second, 
construct the (n+l) x (n+l) matrix with j-th column corresponding to poly- 
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nomial fj,j — 0,...,n, and whose Xij variables are gradually substituted, in 
successive rows, by each respective yij variable. This construction is named 
after Bezout or Dixon and is well-known in the literature, e.g. [CM96], [EM99]. 
A general entry is of the form 

fj{ya(l), ■ ■ ■ ,ya{k-l),ya{k)l, ■ ■ ■,ya{k)t, (11) 
^a{k){t+l)j ■ ■ ■ ■> ^c{k)l„(k)i •^<y{k+l)-i ■ ■ ■ ■> ^a{r)) 

where a :— Tr~^,k = 0,...,r,t = l,...,lk- There is a single first row for 
k — 0, containing all the polynomials in the Xij variables, whereas the last 
row has the same polynomials with all variables substituted by the i/ij. All 
intermediate rows contain the polynomials in a subset of the Xij variables, the 
rest having been substituted by each corresponding yij. The number of rows is 
1 + J2j£[r] = 1 + ^- Lastly, in order to obtain B'^{x, y), we divide the matrix 
determinant by 

r I 

IlIli^^J-y^J)■ (12) 

i=l j=l 

Example 6.11 Let / = (1,2), d = (2,1). If vr = (12), tt' = (21), then m'^ = 
(5, 0),m'^ = (1,1). For both degree vectors, the matrix dimension is 6. To 
obtain B'"{x,y) we construct a 4 x 4 matrix whose j-th column contains 
fj(xi,X2i,X22),fj(yi,X2i,X22), fj(yi, y2i, X22) , /j(yi , Z/21 , Z/22) , fori = 0, ...,3. 
Here Xi (and yi) is a shorthand for xu (and yu)- Then B'^{x,y) contains the 
following monomials in the Xi and variables respectively, 6 in each set of 
variables: 1, xi, 2:21, 2:22, a;ia;2i, a;i,T22, l,yi,y\.yl,y'l.y\. So the final matrix is 
indeed square of dimension 6. More details on this example are provided in 
Section 8. □ 

Lemma 6.12 Let B'^{x,y) = Y.bapx'^y^ where a = {aij),P = {(3ij) G Z", 
i = l,...,r,j = l,...,li. Setai = Ej=i,...,z, "ij; A = J2j=i,...,k Pij, for all a, p. 
Then, < ctj < ,0 < Pi < mj and < ai + Pi < pi, i = 1, . . . ,r. 

Proof. By Lemma 6.10 it suffices to bound ccj. Pi. But ccj is the degree of the 
Xi in the determinant decreased by li in order to account for the division by 
(12). The former equals the product of di with the number of rows where an 
Xij variable appears for any j G [1, /«]. These are the first row, the rows where 
yj are introduced for j G {cr(l), . . . , a{k — 1)} such that a{k) = i, and another 
li — 1 rows when (T{k) — i. The condition on j : 7r(j) < 7r(i) is equivalent 
to r + 1 - 7r(j) > r + 1 - 7r(i), hence a, < + E7r'0)>7r'(i) h = rnj' . 
Similarly, we prove the upper bound on Pi. The rows containing yij for some 
j G [1, li] are those where j G {a(A'; + 1), . . . , o"(r)} : cT{k) = i, another li — 1 
rows when a{k) = i, and the last row. Now, n{j) > k + 1 > k = 7i{i), so 
Pi < -k + di T,nij)>n{i) h = 'ml- Clearly a^, A > 0. □ 
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For generic polynomials, the upper bounds of ai,(3i are attained. The lemma 
thus gives tight bounds on the support of the Bezoutian. 

Theorem 6.13 Assume all defects are zero and B'^{x,y) is defined as above. 
For any tt, {hap) is a square matrix of dimension 

dimXo(m)=f ^ M-.-d^r^^ 
\li,...,lrj (n + 1) 

Furthermore, det{bap) = R{fo, ■ ■ ■ , /n)- 

Proof. First, we show that {ba/s) is square of the desired size. The dimensions 
are given by the number of exponent vectors a, (3 bounded by Lemma 6.12 
which are exactly dimKQ{m'^ ),dimKQ{m'^) respectively. Both rri^jrn" are 
determinantal, hence both of these numbers arc equal to degi?/(n + 1), by 
Theorem 6.7 and Corollary 6.8. R{fo, ■ ■ ■ , fn) divides every nonzero maximal 
minor of the matrix (60/3); cf. [CM96], [EM99, Thm. 3.13]. Since any nonzero 
proper minor has degree < degi?, the determinant of the matrix {bap) is 
nonzero and equals the resultant. □ 

Note that there is not a unique choice of higher differentials in the Weyman 
complexes. We could chase the arrows in a resultant spectral sequence as in 
[GKZ94, Ch. 2, Prop. 5.4] to show that the matrix we propose comes from 
the explicitization of one possible choice. We have followed instead the more 
direct route based on the above property of the Bezoutian in [CM96], which 
uses more elementary tools. 



7 Two examples 

7.1 The bilinear system 

The generic system of 3 bilinear polynomials is 

/o = flO + CilXl + a2X2 + a3XiX2, 
fl = bo + biXi + b2X2 + b3XiX2, 
/2 = Co + CiXi + C2X2 + C3X1X2, 

and has type (1, 1; 1, 1). The degree of the resultant in the coefficients is 
^(^iij = 6. We shall enumerate all 14 possible determinantal formulae in 
the order of decreasing matrix dimension, from 6 to 2, and shall make the 
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corresponding maps explicit. This study goes back to the pioneering work of 
Dixon [Dix08]. 

For TT = (1, 2), Definition 5.2 yields m — (2, 1) and the complex is — > = 

//°(1,0)(?) ^ Ko ^ ^°(2,1) ^ 0. The corresponding determinantal pure 
Sylvester matrix is, when transposed, equal to 

Oo Oi 02 Qa 
bo h &2 &3 

Co Ci C2 C3 

ao 02 di ds 

6o 62 bi 63 

_ Co C2 Ci C3 _ 

with rows corresponding to the input polynomials and the same set multi- 
plied by Xi, whereas the columns arc indexed by 1, Xi, X2, X1X2, xl, x1x2- By 
symmetry, another formula is possible by interchanging the roles of Xi,X2. 
Further formulae are obtained by taking the transpose of these 2 matrices, 
namely with m = (—1,0), where the complex is H'^{—4:, —3) = i?°(2, 1)* — > 
{H^{-S, -2)f = 0)*)^ and with m = (0, -1). Recall the definition of 

duality from equation (2). Sylvester maps, as well as other types of maps, are 
further illustrated in Section 7.2. 

There are additional determinantal Sylvester formulae corresponding to m = 
(2,-1) and m — (—1,2). Their matrices contain the fi and the fi multiplied 

by Xi^ or by a;^^. In the former case, the complex is H^{0, —3)^ = (i7°(0) (8) 
H°{l)*f H\l, -2)3 = {H^{1) H^{0)*f, the transposed matrix is 

0,1 02 03 
bo 61 62 63 

Co Ci C2 C3 

01 03 oq 02 
61 63 bo 62 

_Ci C3 Co C2_ 

and the columns are indexed by 1, xi, X2, ^1X2, x^^, a;r^^2- This construction 
can be verified by hand calculations. 

For m = (1, 1) the complex becomes ^ Ki = i/°(0, 0)('i) © H^{-2, -2) ^ 
Ko — i?°(l, 1) — > 0. The matrix is square, of dimension equal to 4, and 
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hybrid. We compute the two maps by hand; for a larger example see Section 7. 
The foundations for constructing such matrices can be found in [WZ94]. The 
transposed 4x4 determinantal formula is written as follows, by using brackets: 



ao 




02 


O3 


bo 


bi 


b2 


bs 


Co 


Cl 


C2 


C3 



[012] [013] [032] -[123] 



, where [ijk] — det 



Oj Oj Ojfc 

bi bj bk 

Ci Cj Ck 



The matrix rows contain the /j and a rational multiple of the affine toric 
Jacobian, whereas the columns are indexed by 1, Xi, 0:2, a;iX2. This formula is 
obtained in [CDS98] in a more general toric setting. An analogous 4x4 matrix 
corresponds to m = (0, 0). 

There are 4 "partial Bezout" determinantal formulae of dimension 3 x 3 for 
m = (-1,1), (1,-1) and for m = (2,0), (0,2). We omit the details of the 
computation. In the first case, the complex is if^(— 4, — 2) = H^{2,0)* 
H^{—2,0)^ — (i7°(0)* (8) H^{0)y, and a choice of the matrix is, in terms of 
brackets. 



[-02] [-03] + [-12] [-13] 
[0-2] [0-3] + [1-2] [1-3] 
[02-] [12-] + [03-] [13-] 



where [ij—] = det 



tti ttj 



bi bj 



and analogously for the 2x2 brackets [i — k], [—jk]. The columns of this resul- 
tant matrix are indexed by l,xi,xl, which is the support of the 3 Bezoutian 
polynomials filling in the rows. In particular, these polynomials are defined 
for {i,j, k} — {0, 1, 2} in the standard way: 



Bk — det 



fi{xi,X2) Mxi,y2) 

fj{Xl,X2) fj{Xl,y2) 



I {x2-y2)- 



For m = (1, 0) the complex becomes ^ i^i = H'^{-2, -S)^^) = H°{0, 1)* ^ 
Kq — i?°(l,0) — > 0. The corresponding determinantal pure Bezout-type for- 
mula is obtained from the Bezoutian polynomial 



B = det 



Mxi,x2) fo{yi,x2) fo{yi,y2) 
fi{xi,x2) fi(yi,x2) /i(yi,y2) 
f2(xi,x2) f2(yi,x2) /2(yi,y2) 



/ {xi-yi){x2-y2), 
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supported by {1, X2}, {1, yi}- The resultant matrix is given in terms of brackets 
as follows: 

[123] [023] 
-[103] [012] 

7.2 A hybrid determinantal formula 

Assume I — {3,2), d = (2,3). We present explicit formulae which can be 
extrapolated in general, giving an answer to the problem stated in [WZ94, 
p. 578]. We plan to carry this extensively in a future work, but we include 
here the example without proofs as a hint for the interested reader. Our 
Maple program enumerates 30 determinantal vectors m, among which we 
find m^^'^^ = (3, 13) according to Theorem 4.1. 

The minimal matrix dimension is 1320 and is achieved at m = (6, 3) and 
(2, 12). In both cases, P2(m) = 0, whereas Pi(6, 3) = {4} and Pi(2, 12) = {2}. 
This shows that the minimum matrix dimension may occur for some empty 
Pfc, contrary to what one may think. 

Moreover, the degree of the sparse resultant is 6(3^2^32 = 4320. Since 1320 
does not divide 4320, the minimal matrix is not of pure Bezout type; it is 
not of pure Sylvester type either. To specify the cohomologies and the linear 
maps that make the matrix formula explicit we compute, for the degree vector 
m = (6,3) and p = 1, ... ,6 the different values of m — pd: (4,0), (2, —3), 
(0,-6), (-2,-9), (-4,-12), (-6,-15). The complex becomes K2 ^ ^ 
Ki ^ Kq ^ K_i = 0, with nonzero part 

H^{A, 0)(i) e i/^(0, -6)(3) e H^{-Q>, -15)(6) ^ 
^ //°(6,3)(o) e//^(2,-3)(') e//^(-4,-12)(5), 

where we omitted the reference to the space X = P'^ x in the notation of 
the cohomologies. Then dimivTi = 210 + 200 + 910 = 1320 = 840 + 150 + 330 = 
dimiCo- By a slight abuse of notation, let 5a,i3 stand for the restriction of the 
above map to m P"^. Then ^02 = (5o5 = ^5 = by [WZ94, Prop. 2.5] and 
it suffices to study the maps below, of which the first 3 are of pure Sylvester 
type by [WZ94, Prop. 2.6] and the last 3 are of pure Bezout type as those of 
Section 6. These maps can be simplified using the dual cohomologies: 

W{¥\mk - pdk) = H'^-={¥\ {pk -mk)-{n + l- p)dky , 

where p is the critical vector of Definition 2.3. So, we have maps 
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{H%0) H%3y)is) {H°{2) //°(0)*)(') 
//°(2,12)* ^ {H^{0,9)*f 

(//°(0)«)/i'°(3)*)(3) ^//°(6,3) 
//°(2,12)* ^ i/°(6,3) 

//°(2, 12)* ^ (//°(2) (g) H^{Oy)(') 



The resultant matrix (of the previous map in the natural monomial bases) has 
the following aspect, indicated by the row and column dimensions: 





840 


150 


330 










210 


^00 










Soo 








220 




^22 







^20 


^22 





910 




S52 


S55 




-B50 


r>x\ 
-D52 


qT 



where Sij, B^j, B^^ stand for pure Sylvester and Bezout blocks, the latter com- 
ing from a Bezoutian with respect to variables for A; = 1,2, and S^^ repre- 
sents a transposed Sylvester matrix, corresponding to the dual of the Sylvester 
map //0(0,9)6 ^ //0(2,12). 

Let us take a closer look at 522, which denotes both the map and the corre- 
sponding matrix. Let a e N'^, /3 e N'^, be the degree vectors of the elements of 
H^{2), H'^isy respectively, thus |a| < 2, |/?| < 3. Let /, J C {0, . . . , 5}, |7| = 
3, 1 J| =2 express the chosen polynomials according to the cohomology expo- 
nents. Then the entries are given by 



522«®rj,i®^f) 



0, if J ^ /, 

coef(/fe) of<4, if /\ J = {A;}, 



where Tj e if °(0)*, Sj are elements of the respective dual bases of monomials. 
We expect such a construction to be generalizable, but such a proof would be 
part of future work. 

Now take the Bezout maps: The matrix entries are given in (11) for a = (2, 1): 
the entry (i, j), i,j e {0,..., 5} contains /y(xW, . . . , x^^-'\ y^^-'\ y^^)), 
where each x^^^ is a leading subsequence of Xu, X12, Xis, X21, X22', similarly 
with the new variables y^^\ The degree of the determinant, i.e. the Bezoutian, 
is 6, 3, 2, 12 in Xi,X2, yi, y2 respectively and these coefficients fill in the matrix 
i?50. For the Bezout block i?52, consider "partial" Bezoutians defined from 
the 6 polynomials with the exception of those indexed in J, where J, I are as 
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routiiiG 


function 


compjn 


compute the degree vector ni by some specified formula 


a.±xue L V ecs 


enumerate all determinantal formulae 


a.Hsiinis 


compute all possible sums of the Zj's adding to q €z {0, • • • , ^^j— i ^i}' 


coHzero 


test whether H'^i^X^ ui — p^i) vanishes 


coHdim 


compute the dimension of H^i^JC^ ui — 


diniKv 


compute the dimension of Ky i.e. of the corresponding matrix 




iiiiu. ail //t-vecLOlo yicioiiig a puie xDezouu-Lype iuiiiiuia 


f indSyl 


find all m-vectors yielding a pure Sylvester-type formula 


minSyl 


find all m'^-vectors yielding a pure Sylvester-type formula 


hasdeterm 


test whether a determinantal formula exists 



Table 1 

The main functionalities of our software. 



above. Only the xi variables are substituted by new ones, thus yielding a 4 x 4 
matrix. For i?2o, take all polynomials indexed in I and develop the Bezoutian 
with new variables ?/2 from a 3 x 3 matrix. Hence the entries of the Bezout 
blocks have, respectively, degree 6,4,3 in the coefficients of the /«. 



8 Implementation 

We have implemented on Maple V routines for the above operations, in- 
cluding those in Table 1. They are illustrated below and are available in file 
mhomo.mpl through: 

http : //www . di . uoa . gr/~einiris/ index-eng . h tmlj 

Example 3.2 (Cont'd) Recall that I = {l,2),d= (2,3) and let m = (6,3): 

> Ns:=vector([l,2]) : Ds:=vector( [2,3] ) : 

> suinms:=allsums(Ns) : 

> hasdeterm(Ns ,Ds , vector ( [6,3] ) ,summs) ; 

true 

> dimKvCNs ,Ds , vector ( [6,3] ) , summs , 1) ; 

88 

> dimKv(Ns ,Ds , vector ( [6,3] ) , summs, 0) ; 

88 

□ 

Example 5.5 (Cont'd) Recall that / = (2, l,l),rf = (2,2,2), then 5 = 
(1,0,0). The Maple session first computes all 81 pure Sylvester formula by 
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searching the appropriate range of 246 vectors. The smallest formulae are 
shown. 

> Ns : =vector ( [2 , 1 , 1] ) : Ds : =vector ([2,2,2]): 

> ininSyl(Ns,Ds) : 

list of minimal S-matrices: m-vector and Kl,KO-dims 

[[8, 5, 3, 1120, 1080], [8, 3, 5, 1120, 1080], [6, 9, 3, 1200, 1120], 
[6, 3, 9, 1200, 1120], [4, 9, 7, 1440, 1200], [4, 7, 9, 1440, 1200]] 

> allSyl:=findSyl(Ns,Ds) : 

Search of degree vecs from [4,3,3] to [8,9,9] . 
First array [4,7,9]: dimKl=1440, dimK0=1200, 
dimK(-l)=0( should be 0). 

#pure- Sylvester degree vectors= 81 
tried 246, got 81 pure-Sylv formulae [m,dimKl ,dimKO] : 

> sort (convert (%, list) ,sort_fnc) ; 



[[8, 5, 3, 1120, 1080], [8, 3, 5, 1120, 1080], 

[6, 9, 3, 1200, 1120], [6, 3, 9, 1200, 1120], 
[4, 9, 7, 1440, 1200], [4, 7, 9, 1440, 1200], 
[8, 6, 3, 1400, 1260], [8, 3, 6, 1400, 1260], . . .] 

> allDetVecs(Ns,Ds) : 

> allmsrtd := sort (convert (%, list) ,sort_fnc) ; 

From, [-4, -3, -3], to, [9, 10, 10], start at, [-4, -3, -3] 
Tested 1452 m-vectors: assuming Pk's nonempty. 
Found 488 det'l m-vecs, listed with matrix dim: 



allmsrtd := [[6,3,1, 224], [6,1,3, 224], [1,7,5, 224], 
[1,5,7, 224], [3,7,1, 240], [4,7,1, 240], [3,1,7, 240], 
[4,1,7, 240], [6,3,0, 262], [6,0,3, 262], [1,8,5, 262], 
[1,5,8, 262], ...] 

> for i from 1 to nops( allmsrtd ) do print ( 
allmsrtd [i] , Pksets(Ns,Ds, vector ( [allmsrtd [i] [1] , 
allmsrtd [i] [2] , allmsrtd [i] [3] ] ) ) ) : od : 

[6, 3, 1, 224], [[4, 4], [2, 2], [1, 1]] 
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[6, 1, 3, 224], [[4, 4], [1, 1], [2, 2]] 

[1, 7, 5, 224], [[1, 1], [4, 4], [3, 3]] 

[1, 5, 7, 224], [[1, 1], [3, 3], [4, 4]] 

[3, 7, 1, 240], [[2, 2], [4, 4], [1, 1]] 

[4, 7, 1, 240], [[3, 3], [4, 4], [1, 1]] 

[3, 1, 7, 240], [[2, 2], [1, 1], [4, 4]] 

[4, 1, 7, 240], [[3, 3], [1, 1], [4, 4]] 

[6, 3, 0, 262], [[4, 4], [2, 2], [0, NOJNT, 0]] 

The last 2 commands find all 488 determinantal vectors. The smallest formulae 
are indicated (the minimum dimension is 224) and for some we report the 
Pfe's. No determinantal formulae is pure Sylvester. Notice that the assumption 
of empty sets is used only in order to bound the search, but within the 
appropriate range empty P^'s are considered, so no vahd degree vector is 
missed. This is illustrated by the last P3 = marked NOJNT. 

The vectors predicted by Theorem 4.1 are found among those produced above, 
including m^^^s) ^ (6, 5, 3), m^^i^) = (4, 9, 3), m^^^^) ^ (0,9,7). The corre- 
sponding matrix dimensions are 672, 600, and 800. □ 

Example 6.11 (Cont'd) Recall that / = (l,2),(i = (2,1). The only pure 
Bezout formulae are the 2 determinantal formulae of Ex. 6.11, for which we 
have rri^ = (5, 0), m'^' = (1, 1). 

> Ns : =vector ( [1 , 2] ) : Ds : =vector ([2,1]): 

> summs :=allsums(Ns) : 

> f indBez(Ns,Ds,true) ; #not only determinantal 

low — upper hounds, 1st candidate :, [0, 0], [6, 1], [0, 0] 
Searched degree m-vecs for ANY pure Bezout formula. 
Tested 15, found 2 pure-Bezout [m,dimKO,dimKl] : 

{[5, 0, 6, 6], [1, 1, 6, 6]} 
The search examined 15 degree vectors between the shown bounds. It is clear 
that both vectors are determinantal because the matrix dimensions are for 
both 6x6. □ 

9 Further work 

Our results can be generalized to polynomials with scaled supports or with a 
different degree d per polynomial. One question is whether the vectors m', m" 
of Definition 3.3 lead to smaller or larger matrices than m. Notice that certain 
cohomologics, which were nonzero for m, may vanish for m' or m" . We plan to 
complete the description of hybrid determinantal formulae. We would also like 
to answer in general the question stated in §4 of determining a priori the degree 
vectors yielding the smallest determinantal formulae in all possible cases. A 
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problem related to the Sylvester formulae calls for identifying in advance the 
nonzero maximal minor in the matrix, which leads to finding a determinant 
with exact degree in some polynomial. 
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